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METHOD FOR CONTROLLING CONFIGURATION 




OF LASER INDUCED BREAKDOWN AND ABLATION 



Government Rights 

This invention was made with government support 
provided by the Office of Naval Research and the National 
Science Foundation under the terms of No. STC PHY 
8920108. The government has certain rights in the 
invention. 

Field of the Invention 

This invention relates generally to methods 
utilizing lasers for modifying internal and external 
surfaces of material such as by ablation or changing 
properties in structure of materials. This invention may 
be used for a variety of materials. 

Background of the Invention 

Laser induced breakdown of a material causes 
chemical and physical changes, chemical and physical 
breakdown, disintegration, ablation, and vaporization. 
Lasers provide good control for procedures which require 
precision such as inscribing a micro pattern. Pulsed 
rather than continuous beams are more effective for many 
procedures, including medical procedures. A pulsed laser 
beam comprises bursts or pulses of light which are of 
very short duration, for example, on the order of 10 
nanoseconds in duration or less. Typically, these pulses 



are separated by periods of quiescence. The peak power of 
each pulse is relatively high often on the order of 
5 gigawatts and capable of intensity on the order of 10 13 
w/cm 2 . Although the laser beam is focused onto an area 
having a selected diameter, the effect of the beam 
extends beyond the focused area or spot to adversely 
affect peripheral areas adjacent to the spot. Sometimes 

10 the peripheral area affected is several times greater 
than the spot itself. This presents a problem, 
particularly where tissue is affected in a medical 
procedure. In the field of laser machining, current 
lasers using nanosecond pulses cannot produce features 

15 with a high degree of precision and control, particularly 
when nonabsorptive wavelengths are used. 

It is a general object to provide a method to 
localize laser induced breakdown. Another object is to 
provide a method to induce breakdown in a preselected 

20 pattern in a material or on a material. 

Summary of the Invention 

In one aspect the invention provides a method 
for laser induced breakdown of a material with a pulsed 

25 laser beam where the material is characterized by a 

relationship of fluence breakdown threshold (F th ) versus 
laser beam pulse width (T) that exhibits an abrupt, 
rapid, and distinct change or at least a clearly 
detectable and distinct change in slope at a 

30 predetermined laser pulse width value. The method 
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comprises generating a beam of laser pulses in which each 
pulse has a pulse width equal to or less than the 
5 predetermined laser pulse width value. The beam is 
focused to a point at or beneath the surface of a 
material where laser induced breakdown is desired. 

In one aspect, the invention may be understood 
by further defining the predetermined laser pulse width 

10 as follows: the relationship between fluence breakdown 
threshold and laser pulse defines a curve having a first 
portion spanning a range of relatively long (high) pulse 
width where fluence breakdown threshold (F^) varies with 
the square root of pulse width (T* 5 ) . The curve has a 

15 second portion spanning a range of short (low) pulse 

width relative to the first portion. The proportionality 
between fluence breakdown threshold and pulse width 
differ in the first and second portions of the curve and 
the predetermined pulse width is that point along the 

20 curve between its first and second portions. In other 

words, the predetermined pulse width is the point where 
the versus r p relationship no longer applies, and, of 
course, it does not apply for pulse widths shorter than 
the predetermined pulse width. 

25 The scaling of fluence breakdown threshold 

(FtO as a function of pulse width (T) is expressed as F^ 
proportional to the square root of (T* 5 ) is demonstrated 
in the pulse width regime to the nanosecond range. The 
invention provides methods for operating in pulse widths 

30 to the picosecond and femtosecond regime where we have 
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found that the breakdown threshold (Fth) does not vary 
with the square root of pulse width (T^). 
5 Pulse width duration from nanosecond down to 

the femtosecond range is accomplished by generating a 
short optical pulse having a predetermined duration from 
an optical oscillator. Next the short optical pulse is 
stretched in time by a factor of between about 500 and 

10 10,000 to produce a timed stretched optical pulse to be 
amplified. Then, the time stretched optical pulse is 
amplified in a solid state amplifying media. This 
includes combining the time stretched optical pulse with 
an optical pulse generated by a second laser used to pump 

15 the solid state amplifying media. The amplified pulse is 
then recompressed back to its original pulse duration. 

In one embodiment, a laser oscillator generates 
a very short pulse on the order of 10 to 100 femtoseconds 
at a relatively low energy, on the order of 0.001 to 10 

20 nano joules. Then, it is stretched to approximately 100 
picoseconds to 1 nanosecond and 0.001 to 10 nano joules. 
Then, it is amplified to typically on the order of 0.001 
to 1,000 milli joules and 100 picoseconds to 1 nanosecond 
and then recompressed. In its final state it is 10 to 

25 200 femtoseconds and 0.001 to 1,000 milli joules. Although 
the system for generating the pulse may vary, it is 
preferred that the laser medium be sapphire which 
includes a titanium impurity responsible for the lasing 
action. 
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In one aspect, the method of the invention 
provides a laser beam which defines a spot that has a 
5 lateral gaussian profile characterized in that fluence at 
or near the center of the beam spot is greater than the 
threshold fluence whereby the laser induced breakdown is 
ablation of an area within the spot. The maximum 
intensity is at the very center of the beam waist. The 

10 beam waist is the point in the beam where wave-front 

becomes a perfect plane; that is, its radius of curvature 
is infinite. This center is at radius R = o in the x-y 
axis and along the Z axis, Z = 0. This makes it possible 
to damage material in a very small volume Z = 0, R = 0. 

15 Thus it is possible to make features smaller than spot 

size in the x-y focal plane and smaller than the Rayleigh 
range (depth of focus) in the Z axis. It is preferred 
that the pulse width duration be in the femtosecond range 
although pulse duration of higher value may be used so 

20 long as the value is less than the pulse width defined by 
an abrupt or discernable change in slope of fluence 
breakdown threshold versus laser beam pulse width. 

In another aspect, a diaphragm, disk, or mask 
is placed in the path of the beam to block at least a 

25 portion of the beam to cause the beam to assume a desired 
geometric configuration, in still further aspects, 
desired beam configurations are achieved by varying beam 
spot size or through Fourier Transform (FT) pulse shaping 
to cause a special frequency distribution to provide a 

30 geometric shape. 
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It is preferred that the beam have an energy in 
the range of 10 nJ (nano joules) to 1 milli joule and that 
the beam have a fluence in the range of 0.1 J/cm 2 to 100 
J/cm 2 (joules per centimeter square). It is preferred 
that the wavelength be in a range of 200 nm (nanometers) 
to 1 /im (micron) . 

Advantageously, the invention provides a new 
method for determining the optimum pulse width duration 
regime for a specific material and a procedure for using 
such regime to produce a precisely configured cut or void 
in or on a material. For a given material the regime is 
reproducible by the method of the invention. 
Advantageously, very high intensity results from the 
15 method with a modest amount of energy and the spot size 

can be very small. Damage to adjoining area is minimized 
which is particularly important to human and animal 
tissue. 

These and other object features and advantages 
20 of the invention will be become apparent from the 

following description of the preferred embodiments, 
claims, and accompanying drawings. 

Brief Description of the Drawings 

25 Figure 1 is a schematic representation of 

a laser induced breakdown experimental system which 
includes a chirped pulse amplification laser system 
and means for detecting scattered and transmitted 
energy. If the sample is transparent, then transmitted 

30 energy can also be measured. 
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Figure 2 is a plot of scattered energy versus 
incident fluence obtained for an opaque (gold) sample 
5 using the system in Figure 1 operated at 150 femtoseconds 
(fs) pulse duration. 

Figure 3 is a plot of calculated and 
experimental values of threshold fluence versus 
pulse width for gold, with experimental values 

10 obtained for the gold sample using the system of 
Figure 1 operated at 800 nm wavelength. The arrow 
shows the point on the plot where the F^ 
proportional to T s no longer applies, as this 
relationship only holds for pulse widths down to a 

15 certain level as shown by the solid line. 

Figure 4 is a graphical representation of sub- 
spot size ablation/machining in gold based on arbitrary 
units and showing the threshold fluence needed to 
initiate material removal; Rs the spot size of the 

20 incident beam and Ra the radius of the ablated hole in 
the x-y plane. 

Figure 5 is a schematic illustration of a 
beam intensity profile showing that for laser 
micro-machining with ultrafast pulse according to 

25 the invention, only the peak of the beam intensity 
profile exceeds the threshold intensity for ablation/ 
machining. 

Figure 6A and B are schematic illustrations of 
a beam showing the placement of a disk-shaped mask in the 
30 beam path. 
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Figure 7 is a plot of scattered plasma emission 
and transmitted laser pulse as a function of incident 
5 laser pulse energy for a transparent glass sample, Si0 2 . 

Figure 8 is a plot of fluence threshold (F^) 
versus pulse width (T) for the transparent glass sample 
of Figure 7 showing that F^ varying with T h only holds 
for pulse widths down to a certain level as shown 
10 by the solid line. Previous work of others is 
shown in the long pulse width regime (Squares, Smith 
Optical Eng 17, 1978 and Triangles, Stokowski, NBS Spec 
Bui 541, 1978) . 

Figure 9 is a plot of fluence threshold versus 
15 pulse width for corneal tissue, again showing that the 

proportionality between F^ and pulse width follows the T h 
relationship only for pulse widths which are relatively 
long. 

Figures 10 and 11 are plots of plasma 
20 emission versus laser fluence showing that at 170 
(Figure 10) pulse width the F^ is very clearly 
defined compared to 7nm (Figure 11) pulse width where it 
is very unclear. 

Figure 12 is a plot of impact ionization rate 
25 per unit distance determined by experiment and 
theoretical calculation. 

Figures 13A and B are schematic illustrations 
of beam profile along the longitudinal Z axis and 
sharing precise control of damage - dimension along 
30 the Z axis. 



8 



Detailed Descriptio n of the Preferre d RmhorliTiigni-g 
Referring to Figure 1 there is shown an 
5 apparatus for performing tests to determine the laser 

induced breakdown threshold as a function of laser pulse 
width in the nanosecond to femtosecond range using a 
chirped-pulse amplification (CPA) laser system. The basic 
configuration of such a CPA system is described in U.S. 

10 Patent No. 5,235,606 which is assigned to the assignee of 
the present invention and which has inventors in common 
with this present application. U.S. Patent No. 5,235,606 
is incorporated herein by reference in its entirety. 

Chirped-pulse amplification systems have been 

15 described by Jeffrey Squier and Gerard Mourou, two of the 
joint inventors in the present application, in a 
publication entitled Laser Focus World published by 
Pennwell in June of 1992. It is described that CPA 
systems can be roughly divided into four categories. The 

20 first includes the high energy low repetition systems 
such as ND glass lasers with outputs of several joules 
but they may fire less than 1 shot per minute. A second 
category are lasers that have an output of approximately 
1 joule and repetition rates from 1 to 20 hertz. The 

25 third group consists of milli joule level lasers that 
operate at rates ranging from 1 to 10 kilohertz. A 
fourth group of lasers operates at 250 to 350 kilohertz 
and produces a 1 to 2 micro joules per pulse. In 
5,235,606 several solid state amplifying materials are 

30 identified and the invention of 5,235,606 is illustrated 
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using the Alexandrite. The examples below use Ti: Sapphire 
and generally follow the basic process of 5,235,606 with 
5 some variations as described below. 

The illustrative examples described below 
generally pertain to pulse energies less than a 
microjoule and often in the nanojoule range with pulse 
duration in the range of hundreds of picoseconds or less 

10 and the frequency on the order of 1 kilohertz. But these 
examples are merely illustrative and the invention is not 
limited thereby. 

In a basic scheme for CPA, first a short pulse 
is generated. Ideally the pulse from the oscillator is 

15 sufficiently short so that further pulse compression is 
not necessary. After the pulse is produced it is 
stretched by a grating pair arranged to provide positive 
group velocity dispersion. The amount the pulse is 
stretched depends on the amount of amplification. Below 

20 a milli joule, tens of picoseconds are usually sufficient. 
A first stage of amplification typically takes place in 
either a regenerative or a multipass amplifier. In one 
configuration this consists of an optical resonator that 
contains the gain media, a Pockels cell, and a thin film 

25 polarizer. After the regenerative amplification stage 
the pulse can either be recompressed or further 
amplified. The compressor consists of a grating or 
grating pair arranged to provide negative group velocity 
dispersion. Gratings are used in the compressor to 

30 correspond to those in the stretching stage. More 
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particulars of a typical system are described in U.S. 
Patent No. 5,235,606, previously incorporated herein by 
reference. 

An important aspect of the invention is the 
5 development of a characteristic curve of fluence 

breakdown threshold as a function of laser pulse 
width specific to a material. Then identify on such 
curve, the point at which there is an abrupt, or distinct 
and rapid change or at least a discernable change in 
10 slope characteristic of the material. In general it is 
more desirable to operate past this point because of the 
more precise control of the laser induced breakdown (LIB) 
or ablation threshold. 



15 Example 1 - Opaque Material 

Figure 1 shows an experimental setup for 
determining threshold fluence by determining scattered 
energy versus incident fluence and by determining 
threshold fluence versus pulse width. The system 

20 includes means for generating a pulsed laser beam as 
described earlier, and means, typically a lens, for 
collecting emission from the target to a photomultiplier 
tube. Change of transmission through a transparent 
sample is measured with an energy meter. 

25 Figure 2 shows a plot of data obtained from an 

absorbing medium which is gold using 150 fs pulse and 
Figure 3 shows threshold fluence versus pulse width. The 
arrow in Figure 3 identifies the point at which the 
relationship between the threshold fluence and pulse 

30 width varies dramatically. 
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In experimental conditions with wavelength of 
800 nm and 200 fs pulses on gold (Figure 3), the 
5 absorption depth is 275 A with a diffusion length of 50A. 
In the case of nanosecond pulses the diffusion length, 
which is on the order of 10 /im (micron) in diameter, is 
much longer than the absorption depth, resulting in 
thermal diffusion being the limiting factor in feature 

10 size resolution. Empirical evidence for the existence of 
these two regimes is as exhibited in Figure 3, Here both 
experimental and theoretical ablation thresholds are 
plotted as a function of pulse width. An arrow at 
approximately 7 picoseconds pulse width (designated 

15 herein as T or t p ) delineates the point (or region 
closely bounding that point) at which the thermal 
diffusion length (1^) is equal to the absorption depth 
(1/a). It is clear that for a smaller size spot a 
shorter (smaller) pulse is necessary. For spot size on 

20 the order of 1000 A or less, pulse width on the order of 
100 femtoseconds or less will be needed. It is clear 
from the figure that this is the point at which the 
ablation threshold transitions from a slowly varying or 
nearly constant value as a function of pulse width to one 

25 that is dramatically dependent on pulse time. This 

result is surprising. It has been demonstrated that the 
electron thermalization time for laser deposited energy 
in gold is on the order of, or less than, 500 fs and the 
electron-lattice interaction time is 1 ps. The 

30 consequences of this for ultrafast laser pulses is that 
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the energy is contained within the beam spot. In fact 
for energies at or near the threshold for ablation, the 
5 spatial profile of the laser beam will determine the size 
and shape of the region being ablated (Figures 4 and 5). 

Additional experiments were performed to 
measure the amount of recombination light produced as a 
function of the fluence impinging on a gold film. The 

10 technique involved is based upon the experimental setup 
previously described. A basic assumption is that the 
intensity of the light is proportional to the amount of 
material ablated. In Figure 4, the material removed is 
plotted as a function of fluence. A well defined 

15 threshold fluence is observed at which material removal 
is initiated. By having only a small fraction of the 
gaussian beam where the fluence is greater than the 
threshold, the ablated region can be restricted to this 
small area. In Figure 4, R a is the radial position on 

20 the beam where the fluence is at threshold. Ablation, 

then, occurs only within a radius R a . It is evident that 
by properly choosing the incident fluence, the ablated 
spot or hole can in principle be smaller than the spot 
size, Rs. This concept is shown schematically in Figure 

25 5. Although the data for a 150 fs pulse is shown in 

Figure 4, this threshold behavior is exhibited in a wide 
range of pulse widths. However, sub spot size ablation is 
not possible in the longer pulse regimes, due to the 
dominance of thermal diffusion as will be described 

30 below. 
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Additional experiments on opaque materials used 
a 800 nm Ti: Sapphire oscillator whose pulses were 
5 stretched by a grating pair, amplified in a regenerative 
amplifier operating at 1 kHz, and finally recompressed by 
another grating pair. Pulse widths from 7 ns to 100 fs 
were obtained. The beam was focused with a lOx objective, 
implying a theoretical spot size of 3.0 /jm in diameter. A 

10 SEM photo-micrograph of ablated holes obtained in a 

silver film on glass, using a pulse width of 200 fs and a 
pulse energy of 30 nJ (fluence of 0.4 J/cm 2 ) produced two 
holes of diameter approximately 0.3 /xm in diameter. 
Similar results have been obtained in aluminum. 

15 These results suggest that by, producing a 

smaller spot size which is a function of numerical 
aperture and wavelength, even smaller holes can be 
machined. We have demonstrated the ability to generate 
the fourth harmonic (200 nm) using a nonlinear crystal. 

20 Thus by using a stronger objective lens along with the 

200 nm light, holes with diameters of 200 angstroms could 
in principle be formed. 

These examples show that by using femtosecond 
pulses the spatial resolution of the ablation/machining 

25 process can be considerably less than the wavelength of 
the laser radiation used to produce it. The ablated 
holes have an area or diameter less than the area or 
diameter of the spot size. In the special case of 
diffraction limited spot size, the ablated hole has a 

30 size (diameter) less than the fundamental wavelength 
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size. We have produced laser ablated holes with diameters 
less than the spot diameter and with diameters 10% or 
less of the laser beam spot size. For ultrafast pulses 
5 in metals the thermal diffusion length, l^CDt)" 5 (where D 
is the thermal diffusivity and t the pulse time), is 
significantly smaller than the absorption depth (1/a), 
where a is the absorption coefficient for the radiation. 

Those skilled in the art will understand that 

10 the basic method of the invention may be utilized in 
alternative embodiments depending on the desired 
configurations of the induced breakdown. Examples 
include, but are not limited to using a mask in the beam 
path, varying spot size, adjusting focus position by 

15 moving the lens, adjusting laser cavity design, Fourier 
Transform (FT) shaping, using a laser operating mode 
other than TEMoo, and adjusting the Rayleigh range, the 
depth of focus or beam waist. 

The use of a mask is illustrated in Figure 6A 

20 and B. The basic method consists of placing a mask in the 
beam path or on the target itself. If it is desired to 
block a portion of the beam, the mask should be made of 
an opaque material and be suspended in the beam path 
(Figure 6A) alternatively, the mask may be placed on the 

25 target and be absorptive so as to contour the target to 
the shape of the mask (Figure 6B) . 

The varying spot size is accomplished by 
varying the laster f/#, i.e., varying the focal length of 
the lens or input beam size to the lens as by adjustable 

30 diaphragm. 
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Operation in other than the TEMoo mode means 
that higher order transverse modes could be used. This 
5 affects the beam and material as follows: the beam need 
not be circular or gaussian in intensity. The material 
will be ablated corresponding to the beam shape. 

The Rayleigh range (Z axis) may be adjusted by 
varying the beam diameter, where the focal plane is in 
10 the x-y axis. 

Example 2 - Transparent Material 

A series of tests were performed on an Si0 2 
(glass) sample to determine the laser induced breakdown 

15 (LIB) threshold as a function of laser pulse width 

between 150 fs - 7 ns, using a CPA laser system. The 
short pulse laser used was a 10 Hz Ti: Sapphire oscillator 
amplifier system based on the CPA technique. The laser 
pulse was focused by an f = 25 cm lens inside the Si0 2 

20 sample. The Rayleigh length of the focused beam is - 2 
mm. The focused spot size was measured in-situ by a 
microscope objective lens. The measured spot size FWHM 
(full width at half max) was 26 jum in diameter in a 
gaussian mode. The fused silica samples were made from 

25 Corning 7940, with a thickness of 0.15 mm. They were 
optically polished on both sides with a scratch/dig of 
20-10. Each sample was cleaned by methanol before the 
experiment. Thin samples were used in order to avoid the 
complications of self-focusing of the laser pulses in the 

30 bulk. The Si0 2 sample was mounted on a computer 
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controlled motorized X-Y translation stage. Each 
location on the sample was illuminated by the laser only 
5 once . 

Two diagnostics were used to determine the 
breakdown threshold F^. First, the plasma emission from 
the focal region was collected by a lens to a 
photomultiplier tube with appropriate filters. Second, 

10 the change of transmission through the sample was 

measured with an energy meter. (See Figure 1) Visual 
inspection was performed to confirm the breakdown at a 
nanosecond pulse duration. Figure 7 shows typical plasma 
emission and transmitted light signal versus incident 

15 laser energy plots, at a laser pulse width of r p = 300 
fs. It is worth noting that the transmission changed 
slowly at around F^. This can be explained by the 
temporal and spatial behavior of the breakdown with 
ultrashort pulses. Due to the spatial variation of the 

20 intensity, the breakdown will reach threshold at the 
center of the focus, and because of the short pulse 
duration, the generated plasma will stay localized. The 
decrease in transmitted light is due to the reflection, 
scattering, and absorption by the plasma. By assuming a 

25 gaussian profile in both time and space for the laser 

intensity, and further assuming that the avalanche takes 
the entire pulse duration to reach threshold, one can 
show that the transmitted laser energy U t as a function 
of the input energy U is given by 
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U t = kU, U < 

U t = kU^[l + In (U/U^)], U > 

5 

where k is the linear transmission coefficient. The 
solid curve in Figure 7 is plotted using Eq. (1), with 

as a fitting parameter. In contrast, breakdown 
caused by nanosecond laser pulses cuts off the 

10 transmitted beam near the peak of the pulses, indicating 
a different temporal and spatial behavior. 

Figure 8 shows the f luence breakdown threshold 
F^ as a function of laser pulse width. From 7 ns to 
about 10 ps, the breakdown threshold follows the scaling 

15 in the relatively long pulse width regime (triangles and 
squares) are also shown as a comparison - it can be seen 
that the present data is consistent with earlier work 
only in the higher pulse width portion of the curve. When 
the pulse width becomes shorter than a few picoseconds, 

20 the threshold starts to increase. As noted earlier with 
respect to opaque material (metal), this increased 
precision at shorter pulse widths is surprising. A large 
increase in damage threshold accuracy is observed, 
consistent with the multiphoton avalanche breakdown 

25 theory. (See Figures 8 and 9.) It is possible to make 
features smaller than spot size in the x-y focal plane 
and smaller than the Rayleigh range (depth of focus) in 
the longitudinal direction or Z axis. These elements are 
essential to making features smaller than spot size or 

30 Rayleigh range. 
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Example 3 - Tissue 

A series of experiments was performed to 
determine the breakdown threshold of cornea as a function 
of laser pulse width between 150 fs - 7 ns f using a CPA 
5 laser system. As noted earlier, in this CPA laser 

system, laser pulse width can be varied while all other 
experimental parameters (spot size, wavelength, energy, 
etc.) remain unchanged. The laser was focused to a spot 
size (FWHM) of 26 /im in diameter. The plasma emission 

10 was recorded as a function of pulse energy in order to 

determine the tissue damage threshold. Histologic damage 
was also assessed. 

Breakdown thresholds calculated from plasma 
emission data revealed deviations from the scaling law, 

15 a as in the case of metals and glass. As shown in 

Figure 9, the scaling law of the fluence threshold is 
true to about 10 ps, and fails when the pulse shortens to 
less than a few picoseconds. As shown in Figures 10 and 
11, the ablation or LIB threshold varies dramatically at 

20 high (long) pulse width. It is very precise at short 
pulse width. These results were obtained at 770 nm 
wavelengths. The standard deviation of breakdown 
threshold measurements decreased markedly with shorter 
pulses. Analysis also revealed less adjacent 

25 histological damage with pulses less than 10 ps. 

The breakdown threshold for ultrashort pulses 
(< 10 ps) is less than longer pulses and has smaller 
standard deviations. Reduced adjacent histological 
damage to tissue results from the ultrashort laser 

30 pulses. 
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In summary, it has been demonstrated that sub- 
wavelength holes can be machined into metal surfaces 
using femtosecond laser pulses. The effect is physically 
understood in terms of the thermal diffusion length, over 
the time period of the pulse deposition, being less than 
the absorption depth of the incident radiation. The 
interpretation is further based on the hole diameter 
being determined by the lateral gaussian distribution of 
the pulse in relation to the threshold for vaporization 
and ablation. 

Laser induced optical breakdown dielectrics 
consists of three general steps: free electron 
generation and multiplication, plasma heating and 
material deformation or breakdown. Avalanche ionization 
and multiphoton ionization are the two processes 
responsible for the breakdown. The laser induced 
breakdown threshold in dielectric material depends on the 
pulse width of the laser pulses. An empirical scaling 
law of the fluence breakdown threshold as a function of 
the pulse width is given by a 7r p/ or alternatively, 
the intensity breakdown threshold, 1^ = F^/r p . Although 
this scaling law applies in the pulse width regime from 
nanosecond to tens of picoseconds, the invention takes 
advantage of the heretofore unknown regime where 
breakdown threshold does not follow the scaling law when 
suitably short laser pulses are used, such as shorter 
than 7 picoseconds for gold and 10 picoseconds for Si0 2 . 



While not wishing to be held to any particular 
theory, it is thought that the ionization process of a 
solid dielectric illuminated by an intense laser pulse 
5 can be described by the general equation 

dne (t)/dt = i\ (E)ne (t) + (dne (t)/dt) pr - 
(dn e (t)/dt) loss 

10 where n e (t) is the free electron (plasma) density, ti(E) 

is the avalanche coefficient, and E is the electric field 
strength. The second term on the right hand side is the 
photoionization contribution, and the third term is the 
loss due to electron diffusion, recombination, etc. When 

15 the pulse width is in the picosecond regime, the loss of 
the electron is negligible during the duration of the 
short pulse. 

Photoionization contribution can be estimated 
by the tunneling rate. For short pulses, E~ 10 8 V/cm, 

20 the tunneling rate is estimated to be w - 4 x 10 9 sec" 1 , 
which is small compared to that of avalanche, which is 
derived below. However, photoionization can provide the 
initial electrons needed for the avalanche processes at 
short pulse widths. For example, the data shows at 1 ps, 

25 the rms field threshold is about 5 x 10 7 V/cm. The field 
will reach a value of 3.5 x 10 7 V/cm (rms) at 0.5 ps 
before the peak of the pulse, and w - 100 sec" 1 . During a 
At ~ 100 fs period the electron density can reach n e - 
n t [l - exp(-WAt)] ~ 10 11 cm" 3 , where n t ~ 10 22 is the total 

30 initial valence band electron density. 



21 



4 



w 



Neglecting the last two terms there is the case 
of an electron avalanche process, with impact ionization 
5 by primary electrons driven by the laser field. The 
electron density is then given by n e (t) = n 0 x 
exp(n(E)t), where n 0 is the initial free electron 
density. These initial electrons may be generated 
through thermal ionization of shallow traps or 

10 photoionization. When assisted by photoionization at 
short pulse regime, the breakdown is more statistical. 
According to the condition that breakdown occurs when the 
electron density exceeds n^ ~ 10 18 cm' 3 and an initial 
density of n 0 ~ 10 10 cm" 3 , the breakdown condition is then 

15 given by t\t p ~ 18. For the experiment, it is more 
appropriate to use n th ~ 1.6 x 10 2X cm -3 , the plasma 
critical density, hence the threshold is reached when tit p 
~ 30. There is some arbitrariness in the definition of 
plasma density relating to the breakdown threshold. 

20 However, the particular choice of plasma density does not 
change the dependence of threshold as function of pulse 
duration (the scaling law). 

In the experiment, the applied electric field 
is on the order of a few tens of MV/cm and higher. Under 

25 such a high field, the electrons have an average energy 

of - 5 eV, and the electron collision time r is less than 
0.4 fs for electrons with energy U > 5 - 6 eV. Electrons 
will make more than one collision during one period of 
the electric oscillation. Hence the electric field is 

30 essentially a dc field to those high energy electrons. 
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The breakdown field at optical frequencies has been shown 
to correspond to dc breakdown field by the relationship 
E^Cw) = E dc ^(l + w 2 t 2 )\ where w is the optical 
5 frequency and r is the collision time. 

In dc breakdown, the ionization rate per unit 
length, a, is used to describe the avalanche process, 
with ti = a(E)v drift/ where v drlft is the drift velocity of 
electrons. When the electric field is as high as a few 
10 MV/cm, the drift velocity of free electrons is saturated 
and independent of the laser electric field, v^^ ~ 2 x 
10 7 cm/s. 

The ionization rate per unit length of an 
electron is just eE/U L times the probability, P(E), that 

15 the electron has an energy > U i# or a(E) = (eE/UjPfE). 
Denoting E kT , E p, and Ei as threshold fields for electrons 
to overcome the decelerating effects of thermal, phonon, 
and ionization scattering, respectively. Then the 
electric field is negligible, E < E kT , so the 

20 distribution is essentially thermal, P(E) is simply exp(- 
Ui/kT). It has been suggested: P(E) ~ exp( -const/E) for 
E kT < E < E p ; P(E) - exp( -const/E 2 ) at higher fields (E > 
E p ) . Combining the three cases the expression that 
satisfies both low and high field limits: 

25 

a(E) = (eE/Uj exp(-Ei/(E( l+E/E p )+E kT ) . 

This leads to a E 2 r p - l/r p , i.e., the 
fluence threshold will increase for ultrashort laser 
30 pulses when E > /E p Ei is satisfied. 
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Figure 12 is a plot of a as a function of the 
electric field, E. From experimental data, calculated a 
5 according to tit p = 30 and tj = av^^. The solid curve is 
calculated from the above equation, using E A = 30 MV/cm, 
E p = 3.2 MV/cm, and E kT = 0.01 MV/cm. These parameters 
are calculated from U = eEl , where U is the appropriate 
thermal, phonon, and ionization energy, and 1 is the 

10 correspondent energy relation length (l kT = 1 P ~ 5 A, the 
atomic spacing, and 1 ± r 30 A). It shows the same 
saturation as the experimental data. The dashed line is 
corrected by a factor of 1.7, which results in an 
excellent fit with the experimental data. This factor of 

15 1.7 is of relatively minor importance, as it can be due 

to a systematic correction, or because breakdown occurred 
on the surface first, which could have a lower threshold. 
The uncertainty of the saturation value of v^^ also can 
be a factor. The most important aspect is that the shape 

20 (slope) of the curve given by the equation provides 

excellent agreement with the experimental data. Thus, 
the mechanism of laser induced breakdown in fused silica 
(Example 2), using pulses as short as 150 fs and 
wavelength at 780 nm, is likely still dominated by the 

25 avalanche process. 

Opaque and transparent materials have common 
characteristics in the curves of Figures 3, 8 f and 9 each 
begins with F^ versus T H behavior but then distinct 
change from that behavior is evident. From the point of 

30 deviation, each curve is not necessarily the same since 
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the materials differ. The physical characteristics of 
5 each material differ requiring a material specific 
analysis, in the case of sio, (Figure 8) the energy 
deposition mechanism is by dielectric breakdown. The 
optical radiation is releasing electrons by multiphoton 
ionization (M Pi) that are tightly bound and then 
10 accelerating them to higher energies by high field of 
laser^. it is thought that only a small amount of 
relatively high energy electrons exist prior to the laser 
action. The electrons in, turn collide with other bound 
electrons and release them in the avalanching process. In 
15 the case of metal, free electrons are available and 

instantly absorbing and redistributing energy. For any 
material, as the pulses gat shorter, laser induced 
breakdown (LIB) or ablation occurs only in the area where 
the laser intensity exceeds LIB or ablation threshold. 
There is essentially insufficient time for the 
surrounding area to react thermally. As pulses get 
shorter, vapor from the ablated material comes off after 
the deposition of the pulse, rather than during 
deposition, because the pulse duration is so shorty&t 
was also observed that the laser intensity also varies 
along the propagation axis (Figure 13). The beam 
intensity as a function of R and Z expressed as: 



J < £ , E- ) - Wd + 2/3.)* , exp(- 2RVw a .) 
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r 

Y\ u " ' where Z R is the Rayleigh range and is equal to Z P = 7rW, 2 



W 0 is the beam size at the waist (Z = 0). 

We can see that the highest value of the field 
is at Z = R = 0 at the center of the waist. If the 
threshold is precisely defined it is possible to damage 
10 the material precisely at the waist and have a damaged 

volume representing only a fraction of the waist in the R 
direction or in the Z direction. It is very important to 
control precisely the damage threshold or the laser 
intensity fluctuation. 
15 For example, if the damage threshold or the 

laser fluctuation is known within 10% that means that on 
the axis (R = 0) 

20 I(0,Z)/I o = 1/(1 +(Z/Z R ) 2 = .9 

damaged volume can be produced at a distance Z R /3 where 
Z R again is the Rayleigh range. For a beam waist of W Q = 
SQQ?)\^ X then Z R = 7rW„ 2 = nX and the d distance between hole can 



1 

Z R nX. as shown in Figure 13 
7 3 



25 r \ J. 

c 

1 J 
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The maximum intensity is exactly at the center 
of the beam waist (Z = 0, R = 0). For a sharp threshold 
it is possible to damage transparent, dielectric material 
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in a small volume centered around the origin point (Z = 
0, R = 0). The damage would be much smaller than the 
beam waist in the R direction. Small cavities, holes, or 
damage can have dimensions smaller than the Rayleigh 
5 range (Z R ) in the volume of the transparent, dielectric 

material. In another variation, the lens can be moved to 
increase the size of the hole or cavity in the Z 
dimension. In this case, the focal point is essentially 
moved along the Z axis to increase the longitudinal 
10 dimension of the hole or cavity. These features are 
important to the applications described above and to 
related applications such as micro machining, integrated 
circuit manufacture, and encoding data in data storage 
media. 

15 Advantageously, the invention identifies the 

regime where breakdown threshold fluence does not follow 
the scaling law and makes use of such regime to provide 
greater precision of laser induced breakdown, and to 
induce breakdown in a preselected pattern in a material 

20 or on a material. The invention makes it possible to 
operate the laser where the breakdown or ablation 
threshold becomes essentially accurate. The accuracy can 
be clearly seen by the I-bars along the curves of Figures 
8 and 9. The I-bars consistently show lesser deviation 

25 and correspondingly greater accuracy in the regime at or 
below the predetermined pulse width. 

While this invention has been described in 
terms of certain embodiment thereof, it is not intended 
that it be limited to the above description, but rather 

30 only to the extent set forth in the following claims. 
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The embodiments of the invention in which an 
exclusive property or privilege is claimed are defined in 
the appended claims. 
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